Thermodynamic properties of ferromagnetic mixed-spin chain systems 
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Using a combination of high-temperature series expansion, exact diagonalization and quantum 
Monte Carlo, we perform a complementary analysis of the thermodynamic properties of quasi-one- 
dimensional mixed-spin systems with alternating magnetic moments. In addition to explicit series 
expansions for small spin quantum numbers, we present an expansion that allows a direct evaluation 
of the series coefficients as a function of spin quantum numbers. Due to the presence of excitations 
of both acoustic and optical nature, the specific heat of a mixed-spin chain displays a double-peak- 
like structure, which is more pronounced for ferromagnetic than for antiferromagnetic intra-chain 
exchange. We link these results to an analytically solvable half-classical limit. Finally, we extend 
our series expansion to incorporate the single-ion anisotropics relevant for the molecular mixed-spin 
ferromagnetic chain material MnNi(N02)4(ethylenediamine)2, with alternating spins of magnitude 
5/2 and 1. Including a weak inter-chain coupling, we show that the observed susceptibility allows 
for an excellent fit, and the extraction of microscopic exchange parameters. 
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I. INTRODUCTION 

Promoted by the synthesis of various one-dimensional 
(ID) bimetallic molecular magnets, the physics of 
quantum spin chains with mixed magnetic moments 
is of great interest. Typically, quasi- ID mixed- 
spin (MS) compounds display antiferromagnetic (AFM) 
intra-chain exchang e 1 ! 2 ! 3 , which has stimulated theoret- 
ical investigations of AFM MS models using a vari- 
ety of techniques such as spin-wave t\ieory^£&Z&£^, 
variational methods^, the density-matrix renormaliza- 
tion group (DMRG)^, and quantum Monte Carlo 
calculationaS*2*2i2iiS. Interestingly, since a unit cell of the 
MS chain comprises of two different magnetic moments, 
the spectrum will allow for excitation of 'acoustic' as well 
as 'optical' nature^. While not identified unambigu- 
ously in present day experiments, the character of these 
excitations should appear in thermodynamic and other 
observable properties as two independent energy scales*^ 

In addition, and apart from the preceding, what re- 
mains less well studied are MS chains with ferromagnetic 
(FM) intra-chain exchange, which arise in materials of 
recent interest such as MnNi(N02)4(en)2 with en = 
ethylenediamine 13 . This compound is regarded as a 
quasi-lD MS material with spins S = 5/2 and s = 1 
at the Mn and Ni ions, respectively. The susceptibil- 
ity displays an easy-axis anisotropy. At temperatures 
below Tjy = 2.45K, a weak AFM inter-chain coupling 
induces AFM ordering. If this antiferromagnetic order 
is suppressed by a magnetic field of approximately 1.6T, 
the low-temperature specific heat shows a maximum at 
T = 6K and a shoulder at T — 1.5K. These features 
could possibly reflect the two aforementioned character- 
istic energy scales, however for the case of an FM, rather 
than an AFM MS chain. 

Motivated by this, it is the purpose of this paper to 
perform a complementary analysis of the thermodynamic 
properties of FM MS chains, using high-temperature se- 



ries expansion (HTSE), exact diagonalization (ED), and 
quantum Monte Carlo (QMC). In particular, our HTSE 
will be derived for arbitrary alternating spins S and s. 
This will allow not only for a direct comparison with ex- 
perimental data, but also for a study of a gradual tran- 
sit to the 'half-classical' limit S — > oo which is exactly 
solvablo 14 ! 15 ! 16 . Finally, while our ED and QMC data will 
be obtained on systems of the smallest possible mixed- 
spin magnitude, i.e. S — 1 and s = 1/2, this is a case 
also included in the HTSE. In addition, for the sake of 
completeness and to compare with existing literature, we 
will also include results for the AFM case. 

In this paper, we focus on a chain with two kinds 
of spin species, i.e. S and s, arranged alternatingly 
and coupled by a nearest-neighbor Heisenberg exchange. 
Namely, the Hamiltonian reads 
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The subscript i — 1 . . . N refers to the unit cells, and we 
always use periodic boundary conditions. 

In section ITT1 our HTSE approach is detailed. In sec- 
tion II I II we turn to a comparison with QMC and ED 
and the results of the 'half-classical' limit. In section llVl 
we discuss the result of fitting the HTSE to susceptibil- 
ity data obtained for MnNi(N02)4(en)2. Conclusions are 
presented in section Ivl 



II. HIGH-TEMPERATURE SERIES 
EXPANSION 

The HTSE is an expansion in powers of (3J , where 
(3 is the inverse temperature. Here we use the linked- 
clustcr expansion of Ref. yjj In this method, the series 
coefficients for the thermodynamic limit are obtained ex- 
actly from those calculated on finite-size clusters. In gen- 
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eral, this includes the subtraction of contributions from 
a large number of so-called subclusters. In one dimen- 
sion, however, significant simplifications occur due to 
cancellation's*^. This is true also for the MS chain sys- 
tems. That is, in the absence of a magnetic field, the free 
energy F in the thermodynamic limit is represented by 

F/N = Fi(S,a) + F t (a,S) 

-Ft-^S, s) ~ Fi-iia, S) + 0[(PJ) 2 % (2) 

where Ft(S, s) is the free energy of the £-site open-chain 
system described by 

(e-i)/2 

H * = ~ J J2 (S*-s,; + s t -S l+ i) {e:odd) 

i=l 

= He-i — JSi - si (I : even), (3) 

and Ff(s,S) is obtained by exchanging S and s in 
Fi(S, s). Then, a calculation of Tr[(7ig) n ] is needed on a 
finite system, i.e., 

TrH% = '%2(mi...mt\?Q\mi...mt), (4) 

{rrn} 

where m, represents the magnetic quantum numbers at 
site i. We apply He order by order on the ket \mi . . . m e ) . 
This operation yields linear combinations of kets with 
coefficients which are functions of {mi}. To evaluate 
Tr7i 2 ", products of kets of type W™|mi . . . me) are needed 
at most, while in the case of TrH]™ , one can use that 
H e n+1 = H"H" +1 . In order to evaluate this trace, we 
use two different algorithms. 

Method (i) is based on a direct matrix multiplication 
for fixed 5* and s. A linear combination of kets with 
coefficients is regarded as a sparse vector. It is stored as a 
compressed array of non-zero elements and another array 
of their pointers to the kets. These pointers are stored in 
the ascending order so that one can find a needed element 
using binary search in the array. All the operations are 
performed using integers, and thus there is no loss of 
precision. 

Method (ii) is designed for arbitrary spins, which is 
based on an analytic approach to the matrix elements 
in Eq. It has an advantage for large spins because 
the method (i) will fail for very large spins due to time 
and/or memory constraints. After symbolic operations, 
the summation of the form J2m=-s m ™ can ^ e calculated 
analytically for arbitrary n. For example, when n = 2, 
the sum is equal to |s (s + l)(2s + 1). As a result, the 
series coefficients are obtained as an expression valid for 
arbitrary S and s. Naively it may seem that operators of 
type will causes square roots in the matrix elements 
of the Hamiltonian. However, such square roots are ab- 
sent in the final result. In fact, the calculation can be 
carried out disregarding square roots as explained below. 
Introducing the simplified notation, 

\± ni ...±n e ) = (sf) ni . . . (sfr\ mi . . . m e ). (5) 



the initial state is represented by |0. . .0). Spin opera- 
tors act on these states as follows. Suppose | ± n) (n^O) 
represents the state at site i in the notation (JjJJ. Then, 

s z \ ±n) = (m±n)\ ±n), (6) 
s ± \ ±n) = | ±n± 1), (7) 

and 

s T \ ±n) = s T s ± | i«T I) 
= | ± n T 1) x 

{s(s + 1) - (m + n=F l)(m + n)}. (8) 

Note that the norm of \n) is not unity, namely, 

n 

(±n\ ± n) = TT {s(s + 1) - (m + n' =p l)(m + n')} . 

n' = l 

(9) 

Besides the methods (i) and (ii), the contribution to 
the specific heat from the largest cluster is calculated 
separately. Namely, contributions from £-site chain to 
0[(/3J) 2 ^ 2 ] and OpJ) 2 ^ 1 ] have a simple form; with 
notation x = s(s + 1) and X = S(S + 1), for I = 21 it is 
proportional to 21 x l X l and for £ = 21 + 1 to x l X l+1 + 
x l+1 X l . The prefactors of these terms can be determined 
by comparing with those of s = S = 1/2 for any I in 
Ref. ^1 The methods (i) and (ii) are used only for the 
rest of the contribution. 

We have computed the specific heat for the model 
with s = 1/2 and S = 1, up to 29th order using the 
method (i). Furthermore, for arbitrary s and S, the series 
has been calculated up to 11th order using the method 
(ii) and standard symbolic packages^!. By setting S — 
s = 1/2, our model is reduced to the homogeneous spin- 
half Heisenberg chain, and our series agrees with that in 
the literature^. Furthermore, we have also checked that 
the S — ► oo limit of the series agrees with the Taylor 
series of the exact solutionis which will be shown in the 
next section. 



III. LARGE- AND SMALL-SPIN LIMITS 

In this section, we provide evidence for the presence 
of a double-peak-like structure in the specific heat of the 
FM MS chain. We begin by recalling the elementary 
excitations in an MS chain. The dispersion relation of 
the one-magnon excitations in the FM case reads 

=j(s + s± y/S 2 + s 2 + 2Sscos(k)^ , (10) 

and is shown in Fig. ^ f° r 1 ne extreme quantum case, 
s = 1/2 and S = 1. Similar to the AFM casein, the 
spectrum consists of both an acoustic and an optical 
branch reflecting the presence of two different spins in 
a unit cell. These branches indicate two energy scales in 
the thermodynamics of the MS chains. The appearance 
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FIG. 1: The spin-wave dispersion relation for the FM MS 
chain with s = 1/2 and 5=1. The spatial distance of two 
spins of the same species in neighboring unit cells is taken to 
be unity. 

of two one-magnon branches can be qualitatively under- 
stood as follows: At k = n, the excitations correspond to 
an alternating tilting of either the shorter or the longer 
spins as indicated by the arrows in Fig. ^ The magni- 
tude of neighboring spins determines the energies of these 
modes, leading to a splitting of the branches for S ^ s. 
As k is reduced away from k = tt, the magnon eigen- 
states gradually loose these alternating tilting form, and 
near k = become similar in nature to those of uniform 
chains. The low-temperature specific heat as computed 
using linear spin-wave theory is expected to be exacljSS 
to leading order in T. Since uj(k) is proportional to k 2 
at low energies, the specific heat is thus proportional to 
T 1 / 2 for T < J. 

We can establish a connection between the high- 
temperature series for the specific heat and this low- 
temperature scaling using a suitable Pade analysis. In 
extrapolating the high-temperature series for the specific 
heat, information from the ground-state energy, the low- 
energy excitations, and the high-temperature entropy can 
be used as follows^ 3 -: The expansion variable is changed 
from (3 to the internal energy per unit cell, e = {Hint) /N. 
Here, e = for (3 — 0. The power series of e(/3) is thus 
inverted to obtain /3(e). Let S(e) denote the entropy per 
unit cell of the MS chain. From (3 — d6>/de, one obtains 

S = S T=00 + f /3(e')de', (11) 
Jo 

where S T=00 = ln(25 l + 1) + ln(2,s + 1). The low- 
temperature behavior of the specific heat, C oc T 1 / 2 , 
translates into S oc (e— eo) 1 / 3 at e ~ eo for the new series, 
where eo is the ground state energy per site. In Ref. I23L 
only the FM eo is used to extrapolate the specific heat of 
a HTSE for a FM model. Here, we also employ the AFM 
eo, since this additional constraint from the other sign of 
J does not drastically change the final result, but makes 
the extrapolation less sensitive to the used Pade approx- 
imant. We thus obtain the same extrapolation for both 



the FM and the AFM case. Namely, the AFM specific 
heat will be obtained using the substitutions e — > — e 
and (3 J — ► — (3J from the FM case. Then, if the Pade 
approximant in e is applied to 

5 3 (ln(25+l)(2 S +l))" 3 (l -^j (l + -^j , 

the low-temperature behavior and the high-temperature 
entropy are correctly reproduced. For s = 1/2 and 5=1, 
we use the AFM ground-state energy e^ FM = -1.45408 J 
from Ref. 0. The extrapolation is found to be rather in- 
sensitive to errors in e^ FM . For example, an error 10 _4 J 
in e^ FM affects the AFM specific heat only by 4 x 10~ 4 at 
T = 0.1J, and even less at higher temperatures. In the 
following we denote by P[m/n] the rational approximant 
function in e resulting from a polynomial of order m over 
a polynomial of order n. 

In Fig. |2 a comparison is shown between different Pade 
approximants for the low-temperature specific heat for 
the s = 1/2, S = 1 case. We also include results from 
stochastic series expansion QMC simulations^ for both 
the FM and the AFM case for chains of 100 sites. No 
deviations were found to the QMC data for 50 sites, and 
thus we regard these results to represent the thermody- 
namic limit. For the AFM case, the DMRG results from 
Ref. are shown in Fig. [2K as well. In contrast to the 
AFM case, the HTSE result for the low-temperature spe- 
cific heat for the FM case shows large oscillation upon 
increasing the order of the series. Comparison with the 
QMC data shows that the exact solution is located within 
the range of these oscillations. We take the arithmetic 
average of P[12/ll], P[13/12], P[14/13] and P[15/14] as 
the final HTSE result, which in the temperature range 
0.03J < T < 0.45 J has an error of order 5%. 

The specific heat in a larger temperature regime is 
shown in Fig. We find overall good agreement be- 
tween the QMC and the HTSE data, both for the AFM 
and the FM case. In contrast to the homogeneous FM 
spin- 1/2 chain, the FM MS chain displays at least two 
distinct structures in C(T), namely a peak at T ~ 0.54J 
and a shoulder at T ~ 0.25 J. The HTSE result indicates 
the presence of an additional weak shoulder at T ~ 0.1 J, 
which is however difficult to check for using QMC due to 
increasing statistical errors at that low temperatures. 

The features in the intermediate temperature regime 
become more pronounced with increasing order of the se- 
ries expansion, probably because the higher-order poly- 
nomials can reproduce the involved rapid changes more 
accurately. 

We have included in Fig. [3] ED results for small finite 
chains for both the AFM and the FM case. The full spec- 
trum has been calculated for finite chains with up to 7 
unit-cells (14 sites). The ED result for 14 sites agrees 
with the QMC and the HTSE down to temperatures of 
T - 0.5 J for the FM case, and T - 0.3 J for the AFM 
case. Fig. [3] clearly shows that finite-size effects in the 
specific heat are significantly larger in the FM than in 
the AFM case. In the AFM case, the differences between 
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FIG. 2: Low-temperature behavior of the specific heat per 
unit cell of the AFM (a) and FM (b) mixed-spin chain with 
s = 1/2, 5 = 1 using several different Pade approximants of 
the HTSE and from QMC. In (a), DMRG results from Ref. 9 
are also shown. Here, R is the gas constant. 

the 12- and 14-site data are already rather small, with a 
weak shoulder at low temperature signaling the second 
energy scale in C (T) . For the FM case the ED data 
show a strong finite-size shift of the specific heat max- 
imum, which gradually develops into the shoulder near 
T ~ 0.25J, which is found both from the HTSE and 
QMC. 

This suggests that long-range collective modes domi- 
nate the thermodynamics up to higher temperatures for 
the FM than for the AFM case. To shed more light onto 
this difference, let us compare the elementary excitation 
spectra for the FM and the AFM case: (i) At low tem- 
peratures, the specific heat reflects the dispersion of the 
acoustic branch of excitations. In the long- wavelength 
limit, the dispersion in linear spin-wave theory reads 

where the plus (minus) sign represents the FM (AFMiS) 
case. Therefore, the slope of the acoustic branch in 
the AFM case is larger than in the FM case, making 
low-temperature shoulders more pronounced for FM MS 
chains, (ii) The gap between the acoustic and the optical 



FIG. 3: Specific heat per unit cell of the AFM (a) and FM 
(b) mixed-spin chain with s — 1/2, S = 1 as obtained from 
HTSE, QMC and full diagonalizations for small finite chains. 



mode in the FM case is larger than in the AFMLS case, 
so that the splitting of the structures in C(T) is more 
evident in the former case. Apart from (i) and (ii), the 
elementary excitation spectra of the FM and the AFM 
case resemble each other, pointing towards additional ef- 
fects from magnon interactions. 

After analyzing the extreme quantum case with s = 
1/2 and S = 1, we now discuss the specific heat of MS 
chains for larger values of S, keeping s = 1/2 fixed. In 
particular, we want to connect to the exactly solvable^ 
'half-classical' limit, S — > oo. 

Fig. 21 shows the results from the 11th order HTSE 
for s = 1/2 and various values of S. The temperature 
is normalized to |S| J = Jy/ S(S + 1) in order to render 
the 'half-classical' limit finite. In the limit S — > oo the 
specific heat of the chain with FM couplings is identi- 
cal to that of the AFM chain, with a distinct peak at 
T ~ 0.5 |S| J. At T = the specific heat of the 'half- 
classical' model is finite. In the quantum case, however, 
the specific heat vanishes as T — » 0. Using the variable S 
HTSE, we can link both limiting situations as follows: In 
the limit S — » oo, the series coefficients of (j3SJ) n with 
odd n converge to zero. Let us compare the (2n)-th and 
(2n + l)-th order terms at finite S. The ratio of the latter 
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to the former is of 0(t/S), where t = T/(SJ). Hence, if 
the limit S — > oo is taken with t fixed, the latter becomes 
negligible compared to the former as S — > oo. However, 
if S is finite, the two terms are comparable for temper- 
atures below t ~ S^ 1 . Therefore, even if S is large, the 
specific heat deviates from the 'half-classical' behavior 
below t ~ S' -1 , and approaches zero as t — > 0. As a 
result, the specific heat will display a double-peak- like 
structure for large S as well. Summarizing these results 
from the HTSE, there are two different ways of taking 
the large- S and T — > limit, 



lim lim C/R = 1, 

T-+0S-»oo 

lim lim C/R = 0. 



(12) 
(13) 



As seen from Fig. 0] the difference between the FM 
and the AFM cases are most pronounced in the extreme 
quantum limit s = 1/2 and S = 1, and in the low- 
temperature regime. The large- and small-spin limits ex- 
hibit similar structures, suggesting that the specific heat 
of MS chain systems in general show a double-peak-like 
or peak-shoulder structure, both for AFM and FM intra- 
chain exchange, and for any combination of spins. In the 
FM case this structure is more pronounced if \S — s\ is 
large, reflecting the size of the gap between the acoustic 
and the optical mode, similar to the AFM case^ 

We note that the specific heat in the high-temperature 
limit of the AFM case is larger than in the FM case be- 
cause of quantum effects. This is evident from the first 
two terms of the series for the specific heat, which we 
present in Appendix These two terms stem from two- 
site correlations only and dominate the high-temperature 
behavior. Since the total entropy difference between zero 
and infinite temperature is the same in both cases, the 
FM and AFM specific-heat curves therefore have to in- 
tersect at low temperatures. 



IV. FITTING TO EXPERIMENTAL DATA 

We now turn to a comparison to the susceptibility data 
observed on MnNi(N02)4(en)2- In this compound, the 
symmetry around Ni ions is nearly cubic with however a 
fairly large anisotropy at the Mn site to be expected^. 
Hence, we take into account a single-site anisotropy only 
on one of the spins, i.e. S. The (^-factors of the spins S 
and s are represented by G and g, respectively. There- 
fore, the total Hamiltonian reads 



H — Hw 



n 



where 



i=l 

N 

W mag = -h-^(GSj + gsi 



(14) 

(15) 
(16) 
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FIG. 4: The HTSE results for the dependence on S of 
the specific heat per unit cell with s — 1/2 for antiferro- 
magnetic (AFM) and ferromagnetic (FM) coupling, where 
|S| = y/S(S + l). Here, S = oo is from Ref. Q Details 
are in the text. 



Here, we derive the power series of \ m 0D as well as 
/3J up to 0(/3 7 ). When D = and g = G, the series 
coefficients coincide with those in the literatureS^^ as- 
suming a misprinljSi in Ref. l25i Since the contribution 
Hmag is used to evaluate the susceptibility we will only 
consider the case of small Zeeman energies |h| = h — ► 
in the following. The orientation of the magnetic field 
h will be chosen to be either h\\z or h\\x. Most theo- 
retical studies of MS systems are limited to the case of 
D = and g = G, in order to make use of total spin-z, 
i- e - Ylii^i + s f)i conservation. However, for a proper 
comparison to experimental data, D ^ and g ^ G has 
to be accepted, leading to 



[Hint + Wani) Wmag] 7^ 0. 



(17) 



We emphasize, that our HTSE is carried out taking into 
account this non-commutativity. 

Since 3D AFM ordering of MnNi(N0 2 )4(en) 2 below 
Tn = 2.45K signals the presence of a non-negligible inter- 
chain exchange, we will enhance our ID analysis to incor- 
porate this coupling on a phenomenological basis. That 
is to say, fits to the experimental results will be performed 
using an RPA expression 



X 



XlD 



1-JxX 



ID 



(18) 



where xit> is the susceptibility of the pure ID system 
obtained by HTSE and extrapolated by a simple Pade 
approximation (PA). Here, J± effectively models the av- 
erage inter-chain exchange. Figure shows the results of 
our fits of x to experimental data of the susceptibility 28 
with a magnetic field oriented both, perpendicular and 
parallel to the c-axis. Apart from s = 1 and S = 5/2 
we have used g = 2.24 and G = 2 as listed in Ref. [Ij- 
Best fits are obtained for J = 2.8 K, J± = — 0.036 K and 
D = -0.36 K. As estimated from the PA of the HTSE, 
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FIG. 5: Fitting to susceptibility of MnNi(N0 2 )4(en) 2 by 
Eq.QHJ- Here, (/i|A r A/fcB)x is plotted, where /j,b, Na and 
ktB are Bohr magneton, Avogadro's number and Boltzmann 
constant, respectively. 

in the range plotted, the error involved in the theoretical 
curve for xid is within the width of the line. As is ob- 
viously from the figure our theory allows for an excellent 
fit to the experimental data down to T ~4K. Only the 
high-temperature data, for 10K < T < 25K have been 
utilized to set the parameters J, D, and J±. Keeping 
these fixed, the splitting between h||c and h _L c of the 
experimental data tends to become larger than that of 
the theory for T < 4K. 

In Ref.lH values of J - 1.9K and D 0.45K have 

been reported for MnNi(NC>2)4(en)2. These have been 
obtained by fitting a directional average of the suscep- 
tibility with respect to the magnetic field to the 'half- 
classical' limit,— i.e. S — ► oo, applicable for D = 0. 
Moreover the effects of interchain coupling have been ne- 
glected. To compare with this result we may use J± = 
in Eq. 118(1. In fair agreement with Ref. 0, we find 
J ~ 2.4K and D ~ — 0.36K in this case, with a qual- 
ity of the fit however, which is inferior to that shown on 
Fig. [5] 

V. SUMMARY 

We have studied the specific heat and uniform suscep- 
tibility of mixed-spin chain systems using a combination 
of high-temperature series expansion, exact diagonaliza- 
tion, and quantum Monte Carlo techniques. In particu- 
lar, we have contrasted the cases of FM and AFM intra- 



chain exchange. A symbolic high-temperature series has 
been derived for general values of the spin quantum num- 
bers, and including single ion anisotropics. Using this 
series expansion, we were able to extract the microscopic 
model parameters for the quasi-one-dimensional FM 
mixed-spin chain compound MnNi(N02)4(en)2- Com- 
paring our results to the analytically solvable limit S — > 
oo, we found that not only the AFM but also the FM 
case displays a double-peak-like structure in the specific 
heat which is due to the presence of both 'optical' and 
'acoustic' excitations. In fact, we find that for FM intra- 
chain coupling this structure is more pronounced than 
for AFM exchange. The low-temperature specific heat of 
MnNi(N02)4(en)2 shows a weak shoulder around 1.5K if 
magnetic order is suppressed by a finite magnetic field. It 
is thus suggestive to associate this shoulder with the pres- 
ence of 'acoustic' excitations in this compound. However, 
a quantitatively accurate description of the relevant tem- 
perature range for the FM mixed-spin chain with s = 1 
and S — 5/2 in a magnetic field requires further numer- 
ical studies and is left for future investigations. Finally, 
comparing the AFM with the FM case in the extreme 
quantum limit, i.e. s = 1/2, S = 1, we found finite- 
size effects to be more pronounced in the FM than in 
the AFM case. This suggests long-range collective exci- 
tations to be more relevant for the low-temperature ther- 
modynamics of FM mixed-spin chains. 
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APPENDIX A: SERIES DATA 

Using X = S(S + 1) and x = s(s + 1), the specific heat 
series without the single-site anisotropy is given by 

C = Co + C NO (x, X) + C N d(X, x) + Q [(/3J) 12 ] , (Al) 
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2xX 2x 2 X' 
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8xX 643a; 2 X 2 Ax 3 X 
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